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Abstract 



We study dynamics of a scalar field with the non-positively defined Higgs potential 
non-minimally coupled to gravity. In this model the inflationary stage of the Universe 
evolution could change into the quasi-cyclic stage of the Universe evolution with oscil- 



lation behaviour of the Hubble parameter from positive to negative values. Depending 
on the initial conditions the Hubble parameter can perform either one or several cycles 
its evolution. 



1 Introduction 

Models with scalar fields play a central role in current description of the evolution of the 
Universe at the early epoch [TJ El El EH El E] • However, predictions of simplest inflationary 
models are in disagreement with the Planck2013 results [71 E]. The models with the Ricci 
scalar multiplied by a function the scalar field are intensively studied in the inflationary 
cosmology P [HI [TJl [1__J [HI HHl [16]. In the last years, the Higgs inflation has attracted 
a lot of attention pi [131 H EEl [EH] . 

There is also a consideration of theories with nonlocal scalar fields p_7] inspired by the 
string field theory [18] as inflationary models [HI [201 El] • Solutions with a bounce in these 
models have been considered in [221 [23] an d perturbations have been studied in [2H[25]. In 
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one of these models the scalar field is the tachyon of the NSR fermion string and the model 
has the form of a nonlocal Higgs model. Due to the effect of stretching the potential [T9| [20] 
nonlocality destroys the relation between the coupling constant in the potential, the mass 
term and the value of the vacuum energy (cosmological constant) and produces the effective 
Higgs potential with the negative extra cosmological term [26J. This effective Higgs potential 
is non-positively defined. It turns out that in the case of a minimal coupling non-positivity 
of a potential leads to the change of modes: from expansion to contraction [26]. In the cyclic 
model [27] the ekpyrotic phase is generated by a scalar field with a negative potential and 
one has a bounce from a phase of contraction to a phase of expansion, see [28] for discussion 
of the Planck2013 results support of this model. 

In this paper, we consider the induced gravity model, which action includes the term 
proportional to the Ricci scalar multiplied by the square of the scalar field. The model 
incorporates only one scalar field, in other words, the scalar field which couples to the Ricci 
scalar is the same field whose potential drives inflation. We consider the non-positively 
defined potential, which is the Higgs potential minus a positive constant. The characteristic 
property of models with non-positively defined potential is the existence of the forbidden 
domain of the phase plane, which corresponds to non-real Hubble parameter. By numerical 
calculations we explore the case, when the forbidden domain is separated on the two parts and 
show by numerical calculations that the phase trajectories are attracting to the boundary 
of the forbidden domain, touch it and go to infinity. During this evolution the Hubble 
parameter oscillates from positive to negative values. This quasi-cyclic behaviors end up 
on a point of the boundary of the forbidden region and after this the evolution changes 
and completes by monotonically contracting stage. The phenomenon of a quasi-cyclic stage 
of the Universe evolution is a subject of our paper. Note that we got some results for an 
arbitrary non-positively defined potential as well. 

The plan is the following. We start, Section 2, from a short remind of the action and 
equations for the induced gravity cosmological models. We consider dynamics in the models 
with an non-positively defined Higgs potential. We show that there are solutions that wind 
up on the forbidden region on the phase plane. In Section 3, we discuss the model in the 
Jordan and Einstein frames. In Section 4, we study the special case £ = —1/6, in which 
some equations, considering in Section 2, have no sense. Finally Section 5 is devoted to the 
conclusion. 



2 Classical dynamics in cosmological models with non- 
minimal coupling and non-positively defined poten- 



Cosmological models in framework of a non-minimally coupled theory are actively studied [91 
HOI HH EH EOl EH E21 ESI [Ml ES] (see also [3S1 [37] and references therein). Let us consider 
the the induced gravity models described by the following action [TT], [33] : 



tials 




(i) 
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where £ is the no n- minimal coupling constant. 

We study the evolution of a homogeneous scalar field on a spatially flat Friedmann- 
Robert son- Walker (FRW) universe with the metric 

ds 2 = -dt 2 + a 2 it) [dx\ + dx\ + dxQ . 

The Einstein equations in the FRW metric look as follows: 



H* = L£>+±IZ\ - 2H Z (2 ) 




3/7* + 2/7 = - 2^ - 4^ - ^ ^ j + ^(«, (3) 

where differentiation with respect to time i is denoted by a dot, and the Hubble parameter 
is the logarithmic derivative of the scale factor: H = a/a. 

The equation for the field got by the variation of the action over the field is 

4> + 3H0 + V - 6£0 (2H 2 + = 0. (4) 

Combining equations (J2J) (J3J) , we obtain the following system of the first order differential 
equations 

= V>, (5) 

t = - 3^ - j + jrrm m<l)) ~ (6) 

fr = *h± , YM _ J^Lh 2 - 1 + 2e (±\ 2 (7) 

(1 + 600 (1 + 600 l + 6£ 2£(1 + 6OV0/' 

We assume that £ 7^ —1/6. The case £ = —1/6 will be considered separately in Section HJ 
All numerical calculations have been made for £ > 0. 
Equation (J2J) is the quadratic equation for H: 

T 2 , 




# 2 + 2H- = 0, (8) 

3£0 2 6£ 1 f 



and has the following solutions 



The function if is a continuous function, so, if V(0) > for all 0, then evolution of 
the Universe in such model is described either only H or only H + . It depends on initial 
conditions. If V(0) is not positively defined, then it is possible that a part of evolution is 
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described by H_, whereas another part — by H + . The simplest way to get a non-positively 
defined potential from the known positively defined one is to subtract a positive constant. 

Let us calculate a derivative of the expression H + ^ over time. Using system (jSJ)-©, it 
is easy to show that 



db 



2f(6£ + l)0 : 



8£(V - 3£4> 2 H 2 ) - (1 + 6£ + 24^ 2 )^ 2 - 36£ 



(10) 

We can eliminate the potential V from this formula, using equation (121) . Also, it is useful 
to write the Hubble parameter as 

H= - ^- + Q. 



Now equation (11 01) has the following form 

On the (<f), cb) plane there is the boundary of the forbidden region. At any point of this 
boundary Q = 0, so we get that always Q ^ and Q < at ib ^ 0. This means that if a 
solution for system dSD - © with H + touches this boundary at ib ^ 0, then it change to 
solutions with H_ can not touch this boundary at ib ^ 0. 

From equations ([2]) and ([3]) we get the following equation, which form does not depend 
on the potential: 

"*-^H- 3 T i (?)'- (12) 

Substituting the value of H on the boundary of the forbidden region, we get the following 
equation 

It means the boundary of the forbidden region is not a solution for system (jHJ) - © at any 
values of parameters. Only trivial solutions with ib — can belong to this boundary. The 
above-mentioned results are valid for any non-positively defined potential. 

Models of this type in the case of the positively defined Higgs potential have been con- 
sidered in a lot of works (see, for example, [121 E2 ttH EE]). I n this paper, dynamics in the 
model with the non-positively defined Higgs potential 

V H (<b) = e -{<P 2 -b 2 ) 2 -A, A>0, (14) 

where b and A are constants, is studied. 

In the case of the V#(0) potential, the boundary of the forbidden region (Q = 0) has the 
following form 

(l + 6O0 2 = 2A-|(0 2 -6 2 ) 2 . (15) 
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Most attention will be paid to the case when A < eb A /A. In this case, the forbidden region 
on the (4>, (j>) plane consists of two separate parts. 

Numerical calculations give the following phase trajectory for system As we see 

there are two stages of evolution on the phase plane. The first stage corresponds to some 
kind of the cyclic Universe, then the phase trajectory reaches the boundary of the forbidden 
region and the second stage starts, during this stage the Hubble parameter rapidly decreases 
and the Universe contracts (Seed]). 



•0 




Figure 1: The solution of system ©-© at A = 0.05. We choose b = 1, e = 10, f = 10. 
The initial conditions are 4>o = 2, ip = 0, H is calculated by (EJ) with sing "+". On 
the left picture, the phase diagram is present. It is magenta at the corresponding value of 
H = H + > 0, blue at H = H + < 0, red at H = H_ > 0, and green at H = H_ < brown. 
On this picture, brown dashed line corresponds to H + = 0, brown dashed line with long 
dashes corresponds to if_ = 0, black line is the boundary of the forbidden region. On the 
right picture, the Hubble parameter as function of the cosmic time is present. Brown color 
means that H = H + , whereas H = H-. is drown in dark green color. 

The fact which gives the change of the stages is that when the phase trajectory reaches the 
boundary of the forbidden region the square root in the equation for the Hubble parameter 
gets the "minus" sign. 

Note that the obtained behaviors of the solutions is essentially different from the behaviors 
of the solutions in the case A = 0. At A = the forbidden region is absent and the point 
(H = 0, cf) = b, ip = 0) is an attractive fixed point and the Hubble parameter is positive at 
any point (see Fig |5J). 

The cyclic stage can be interpreted as a heating stage. To estimate the period of this 
reheating we can estimate a period of oscillations for a limit trajectory - a trajectory which 
coincides with a boundary of the forbidden region (1 + 6£)^- = — V ((/)), we remind that this 
trajectory is not a solution for the system of equations under consideration. A solution to 
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Figure 2: The solution of system (JH])-© at A = 0. We choose b = 1, e = 10, £ = 10. The 
initial conditions are O — 1-5, ipo — 0, H is calculated by (jHJ) with sing "+". The Hubble 
parameter is always H + . 



equation = — j^V'(<f)) with negative energy E = A' — where A' = has the form 
(see, for example, [38] ) 

= Adn(ttt + c,k) (16) 

where dn(«, k) is the elliptic function of argument u and modulus k, A is the amplitude and 
c is the phase. The frequency Q and the modulus k of the elliptic function are obtained from 
parameters e', b of the equation and depend on the amplitude A: 



A 1 = 6M 1 + a 1 + 



Q 2 EE A 



4£ 
i 7 ^ 



EE 2 1- 



Parameter c is defined by the relation 

0(0) = Adn(c, fc). 



(17) 

(18) 
(19) 

(20) 



Now we estimate a period of function <p(t) (it will be of the same order as the period for 
the function H). Then we will also have an estimation for a period of H. Let us suppose 
the period has the standard form 
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and substitute ft in the form (ITBl . We get 



T = 27tv / 2W il 6 ^ ■ (22) 
2^ + 6% 



Let us estimate the value of the period. For this purpose let us note that a non-minimal 
coupling turns the Planck mass into a dynamical quantity. The Planck mass is related to 
the vacuum expectation value of the potential b as follows [TT] : 

§ = > 2 < 23 > 

This gives us an estimation for £ 

Z = jr4- (24) 

To estimate the value of the constant A let us note that we are considering the case when 
the forbidden area consists of two separate parts. It means that A < let us suppose that 
A ~ £j£ From A ~ it follows that 



e& 2 

eA = — . (25) 



Then for a period it is the following estimation 

3M„ 2 , 

T = -jg. (26) 

For b = 240GeV and ma = y/sb = 126GeV [12], [IS] we get £ ~ 1 and for a period we have 
an estimation 

T ~ lO^GeV -1 ~ 1Q- U s. (27) 



3 Action and equation of motion in the Einstein frame 

Once a modified gravity theory is recast into its scalar-tensor presentation, it immediately 
follows that both the Jordan frame (where the scalar fields non-minimally couples to gravity) 
and the Einstein one (which has the same form as that of Einstein gravity with minimally 
coupled scalar fields) are available. These two frames are related by conformal transformation 

^ = tffl{? 9 = ^9 {E \ (28) 

where we denote the metric in the Jordan frame by g^, while the one in the Einstein frame 
is labeled as g^v . In the following, we denote all quantities in the Einstein frame by adding 
a tag ^ to the corresponding ones in the other frame, in order to avoid confusion. 

Let us consider the conformal transformation ( 128]) . Using g^ v = ft" obtains 
the relationship between the Ricci scalars in the two frames [36J 

R = ft" 2 [R (E) - 6 (D^ In ft + g^ E ^P In ftVf } In ft)] . (29) 
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Inserting this into action ([T]), one identifies the conformal factor as 

C 



n 



(30) 



where C = 1/(«VD> where k 2 = 8tt/M p1 2 , the Planck mass being M P1 ~ 1.22 x 10 19 GeV. 
This results to the following action in the Einstein frame: 



S 



d A xy/-gW 



Ik 2 



R (E) _ 12 + C 2 au(E) d 

K 20 2 9 » 



(31) 



Note that the term proportional to ln(0) has been discarded as a total derivative. One 
can derive the equations of motion in the Einstein frame by varying this action. 

As is well known [36], conformally flat metrics are mapped into each other. The FRW 
metric is conformally flat, so starting from the FRW metric in the Jordan frame we obtain 
the corresponding FRW metric in the Einstein one. This leads us to directly start from an 
FRW metric with cosmic time in the Einstein frame 



dt 2 E + b E {t E )8 i jdx l dx 1 * 



(32) 



where in d&E and o>e^Pe)i the index E denotes the corresponding quantities in the Einstein 
frame. We get 



d log a E 
dtE 



n h 



n 



dt E = gdt, 



C <H + * 



a E = . 



(33) 



±- 



C 



A 



1 + 



The Einstein frame is convenient for calculating the spectral index n s for primordial 
perturbations [UJ . It was shown that the equations (J2J) and ([6]) in this frame can be written 

as 



H 2 E 



d 2 f + 3Rh d( P 
dt E dt e 



K 

~3~ 



U(<p) + 



dip 
dtE 



(35) 
(36) 



where the scalar field (p and the potential U(<p) are 

dip 



l + 6£ 
& 2 ' 



U(ip) 



v{4>) 



(37) 



In slow-roll approximation the Einstein frame equations are 

.2 



I 2 



3H, 



dip 
dtE 



-U'(ip). 
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In other words, it has been assumed that the standard Einstein-frame "slow-roll" approxi- 
mations may be made. This approximate equations in the Jordan frame look as 

2 




\ [4^(0) - . 



Yet the assumption, U((p) 2> ^(dip/dtE ) 2 , is equivalent to V((p) 3> |</> 2 (1 + 6£), and thus it 
remains consistent further to simplify the field equations during slow-roll as: 



" -'<o 2 ' 



3H * - (l + 6O0 [4l/(0) " 0V " (0)] ' (39) 
Using this approximation it is possible to obtain the values for the slow-parameters and the 
spectral index n s as it has been done in pi)] . In the case of non-positively defined potential 
the results will be qualitative the same. One can obtain both "chaotic" inflation and "new" 
inflation. 



4 The special case £ = —1/6 

This case is of a special interest due to the fact that if £ = —1/6 then equations ([5j) (JTJ) have 
no sense. To consider this case we put £ = — | into equations (J2])-(jl]). These equations get 
the following form 




3H> + 2H=\Zj -2r- iH r--V W , (41) 

<j> + 3H<j) + V' + <p (2H 2 + = 0. (42) 
Combining equations (jiOj) and fj4T|) we get 




H = 2\^\ +2H^ + (43) 



After, from equation (|42p we get the condition 

4V((f)) - <f>V'(<f>) = 0. (44) 

This means that a solution with a nontrivial 0(t) exists if and only if the potential 
V{4>) = A<j) A . Constant A must be non-positive to satisfy the requirement that H is real. 
Equation (|4"3]) is an identity in this case. We get that for any function there is function H 
and they are solutions of equations. 
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5 Conclusion 



In this paper, we have considered the with non-positively defined potential, which is the 
Higgs potential minus a positive constant. The characteristic property of models with non- 
positively defined potential is the existence of the forbidden domain of the phase plane, 
which corresponds to non-real Hubble parameter. 

We have shown that behaviors of the Hubble parameter in with non-positively defined 
Higgs potential and Higgs potential are essentially different. For the Higgs potential we have 
attractive fixed points <fi = ±a, which corresponds to minimums of the potentials. If the 
subtract a positive constant from the Higgs potential, then these fixed points appear in the 
forbidden domains. In the case of the non-positively defined Higgs potential the boundary 
of the forbidden region consists of two closed curves (if there is a special relation between 
a cosmological constant and parameters of the potential). We have explored this case and 
shown numerically that the phase trajectories are attracting to the boundary of the forbidden 
domain, touch it and go to infinity. 

In this paper, we have considered the induced gravity model, which are very popular 
in cosmology [HI EQl EE]- Note that cosmological models with R, multiplied by £</> 2 + 
J, where J is a constant, are considered as a real candidate to describe the inflationary 
scenario p21 [13j [15j [16], which is fully consistent with the Planck constraints [7]. We plan 
to consider this model in future investigations. 
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